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Abstract 

Probability distributions which can be obtained from superpositions of Gaussian distributions 



<. 

£ — . . of different variances v = a 2 play a favored role in quantum theory and financial markets. Such 

<n : 

superpositions need not necessarily obey the Chapman-Kolmogorov semigroup relation for Marko- 
Th" 1 ' 

' vian processes because they may introduce memory effects. We derive the general form of the 

smearing distributions in v which do not destroy the semigroup property. The smearing technique 
has two immediate applications. It permits simplifying the system of Kramers-Moyal equations 
for smeared and unsmeared conditional probabilities, and can be conveniently implemented in the 
path integral calculus. In many cases, the superposition of path integrals can be evaluated much 

o: 

easier than the initial path integral. Three simple examples are presented, and it is shown how the 
technique is extended to quantum mechanics. 
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I. INTRODUCTION 



Path integrals are a powerful tool in diverse areas of physics, both computationally and 
conceptually. They often provide the easiest route to derivation of perturbative expansions 
as well as an excellent framework for nonperturbative analysis [1] . One of the key properties 
of path integrals in statistical physics is that the related time evolution of the conditional 
probabilities P(x b ,t b \x a ,t a ) fulfills the Chapman-Kolmogorov (C-K) equation for continuous 
Markovian processes 



P(x b ,t b \x a ,t a ) 



dx P(x b , t b \x, t)P(x, t\x a , t a ) 



(1) 



Conversely, any probability satisfying this equation possesses a path integral representation, 

uation ([!]) also serves as a basis for deriving a 



as shown by Kac and Feynman [2j, 
Fokker-Planck time evolution equation 



i. 



, |5|, [6] for P(x b , t b \x a , t a ) from either the stochastic 
differential equation obeyed by the variable x or the Hamiltonian driving the time evolution 
of P(x b ,t b \x a ,t a ). Such equations are used to explain many different physical phenomena, 
for example turbulence 7| or epitaxial growth 8j. In information theory they serve as a 
tool for modeling various queueing processes [9(, while in mathematical finance they are 



conveniently applied in theory of option pricing for efficient markets 



an 



ii 



12|. 



A trivial property of P(x b , t b \x a , t a ) satisfying the C-K equation (CQ) is the initial condition 



P{Xbi ta | X a , t a ) S(x b X a ). 



(2) 



The right-hand side can be written as a scalar product of Dirac's bra and ket states (x b \ and 
\x a ) as 

P{x bl t a \x ai ta) (x^|x a ). (3) 

In many practical applications one encounters conditional probabilities formulated as a 
superposition of path integrals in which the Hamiltonians H are rescaled by a factor v, i.e., 



P(x b , t b \x a , t a ) 



> i>x(t b )=x b 

dv u(v, t ba ) j VxVp e 

J x(t a )=X a 



fta dr[ipx-vH(p,x)} 



(4) 



where u(v,t ba ) is some positive, continuous and normalizable smearing function of v > 
and t ba = t b — t a > 0. For instance, probability distributions which can be obtained from 
superpositions of Gaussian distributions of different volatilities v = a 2 play an important 
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role in financial markets [l[ [h]]. Such smearing distributions show up also in nonpertur- 
bative approximations to quantum statistical partition functions [131], in systems with time 
reparametrization invariance Q, Q, in polymer physics Q, 3], in superstatistics [if]], etc. 
Whenever smeared path integrals fulfill the C-K equation, the Feynman-Kac formula ensures 
that such superpositions can themselves be written as a path integral without smearing 



P(x b ,t b \x a ,t a ) = / """ VxVpe^ dT l ip± - S ^], (5) 

Jx(t„)=x„ 



rx(tb)=Xb 

lx{t a )=X a 

with a new Hamiltonian H (p, x) given by 

e~ # drR{p > x) = / dv u(v, t ba )e- # d ™ H ^ . (6) 



In general, the smeared path integral (jl]) does not conserve the C-K equation. Physically 
this implies that the superposition (j4j) introduces memory into the system. The aim of 
this note is to find conditions for the smearing distributions where this is avoided. Various 
physical consequences will be derived from this. 

The paper is organized as follows. Section [III starts with a warm up example of smearing 
distributions for a Gaussian conditional probability. In Section HTT1 we derive the most general 
class of continuous smearing distributions fulfilling C-K equation. Section [TV] is devoted to 
a construction of the Hamiltonian H and to a discussion of related issues. In Section [V] we 
show how the outlined path integral representation can be physically interpreted in terms 
of two coupled stochastic processes. Section [VI] discusses three specific smeared systems 
without memory. In particular we show how the explicit knowledge of H (p, x) can streamline 
practical calculations of numerous path integrals. Various remarks and generalizations are 
proposed in the concluding Section rVHl For reader's convenience we include two appendices 
where we perform some finer mathematical manipulations needed in Section [V[ 

II. SMEARING OF A GAUSSIAN DISTRIBUTION 

Our goal is to find the most general form of uj(v, t) fulfilling the C-K relation ([1]). Let us 
first illustrate what we want to achieve by smearing out a simple Gaussian system whose 
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Hamiltonian is H = vp 2 /2 leading to a conditional probability 

r-x(t b )=x b 



P v (x b} t b \xa,t a ) = [ X " Xb VxVpe^ dT l ip± - vp2/2 ] 

J X{t a ) = X a 



' e -(Xb-%a) 2 /{2vt ba ) ^ 



\J1lXVt ha 

where t ba = t b — t a . This obeys the Fokker-Planck equation 

dt b Pv(x b ,t b \x a ,t a ) = — d Xb P v (x b ,t b \x a ,t a ). (8) 

which can be solved explicitly, including the initial condition ([2]) with the help of the differ- 
ential operator for the momentum p = —id x as 

P v (x b ,t b \x a} t a ) = e- tb ^ 2/2 5{x b -x a ), (9) 

or, using the Dirac bra and ket states in ([3]), as 

P v (x b ,t b \x a ,t a ) = {x b \e~ tb ^ 2 ' 2 \x a ). (10) 

Due to the completeness relation J dx\x)(x\ = 1, this expression obviously satisfies the C-K 
equation (CQ). 

Let us now assume that cu(v, t) can be written as a Fourier transform of the form [compare 
also Ref. Q]: 

w(v,t te )= r ^e^- H ^/^ } (11) 

where H u (£) can be viewed as the Hamiltonian affiliated with distribution u. Then the 
smeared transition probability has the integral representation 



P{x b ,t b \x a ,t a ) = 7T--\ x b 



e -E u {£/t ba )t ba 



f - haP 2 /2 



X a ). (12) 



Assuming that H u (£) is a regular function which becomes infinite on a large semicircle in 
the upper half plane we can use the residue theorem to evaluate this as 

P(x b ,t b \x a ,t a ) = (x b \e- H ^ 2 ^\x a ) . (13) 



This can be written path integral 

P(x b ,t b \x a ,t a ) = r (tb) ~ Xb VxVpe^ dT l ip± - H ^ 2 ^]. (14) 

J x(t a )=X a 
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Thus the Hamiltonian H(p, x) of the smeared system in Eq. ([5]) is simply equal to H UJ (p 2 /2). 

In the following we shall generalize this treatment to non-Gaussian Hamiltonians in 
Eq. ©. 



III. SMEARING OF A GENERAL DISTRIBUTION 

We now embark on finding the most general smearing function u(v, t ao ) to guarantee the 
C-K relation for the smeared expression (jlj). Replacing the probabilities in jT]) by (H]), we 
bring the right-hand side to the explicit form 

OO POO POO f%(tb) =x b t 

dv'u(v',t') / dv" u(v",t") dx VxVpe^ dT{ip± - v " H) 

JO J -oo Jx(t c )=x 



% (^a ) — %a 

where t' = t ca = t c — t a and t" = tj, c = t — t c . In the second line we have used the 
substitution t"v" + t'v' = tv and t"v" — t'v' = (. Comparing the right-hand side of ( [151) with 
the left-hand side of (jlj) expressed in the smeared form (jlj), we obtain an integral equation 
for the smearing function uj(v,t): 



tc 



-tv 



( tv -0 .A f( tv + A 2t't"u{v,t) 



It' \ It" 1 / t 



l^—^.,q u ^—^,t" = (16) 



Setting z = tv/t" and going over to an integration variable z' = {tv + ()/2t", this becomes 



dz'u(z\t")u(j(z-z'),?) = , (17) 



dz' u(z\t) auola(z — z'),-\ = bu(bz, ^-), (18) 



or equivalently 



a) \ bj 

with positive real a, b satisfying 1 + 1/a = 1/b. Since the left-hand side is a convolution 
integral, the solution of this equation is found by a Laplace transformation. Defining 



o 



poo 

(£,*) = / dze-e'u(z,t), Re£>0, (19) 
Jo 



we may reduce Eq. (fTBj) to the functional equation 

^•'Ki-i) =<!•!)■ (20) 

It should be stressed that due to assumed normalizability and positivity of u(v,t), the 
smearing distribution is always Laplace transformable. The substitution a = 1/a transforms 
flU} to 

u)(£, t) u(a£, at) = <2(£ + a£,t + at) . (21) 

By considering only real £, Eq. (121]) can be solved by method of iterations familiar from the 
theory of functional equations [18{]. Assume for a moment that a is a positive integer, say 
n, then successive iterations of Eq. (j2Tj) give 

Co(n^nt) = [^,t)] n - (22) 

Now let r = m/n be a positive rational number (m and n positive integers) and ( and r 
arbitrary positive real numbers. Then, for £ = r( = (m/n)( and t = rr = (jn/n)r, we have 
n£ = m( and nt = mr, so that Eq. (|22p yields 

u)(£,t) = cD(rCrr) = [^(C,r)] r , (23) 

for all positive ( and r and all positive rationals r. Assuming that u is continuous, we may 
extend the Eq. ff23l to all positive real r. It is then solved by 

w(f,t) = u(t ■ £/t,t ■ 1) 

= [u,(£/t, l)] 1 = m/t)}\ t>0, (24) 

where G(x) is any continuous function of x. The above derivation is meaningless for t = 0. 
In this case we must instead of f[2~4"|) consider 

<2(£,0) = c^-U-0) = [u>(l,0)]t = k^. (25) 

The constant k is determined by the initial value of the smearing distribution u(v,t). Thus 
Eq. (|25|) implies that 

lim^a;^,^) = 6{v + log + log k) 

= 5 + (v + logK) . (26) 
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Note that Eq. (1241) implies positivity of w(£,£) for all t and £, hence G(x) must also be 
positive for all x. This allows us to write 

[G(t/t)Y = e~ Fmt > (27) 

where F(x) is some continuous function of x. The associated inverse Laplace transform gives 
then the complete solution for u(v,t). Usually, the inverse Laplace transform is expressed 



complex Bromwich integral [19|] [compare (TTTTi ] : 

co( V ,t) = — de^wce.o, (28) 

where the real constant 7 is such that it exceeds the real part of all the singularities of 
£>(£,£). For our purpose it will be preferable to use, instead of (1281) . the inversion formula 
due to Post |2o|: 

k /U\ k+1 



"(•'1 = &t(;J WW.- < 29 > 

For practical calculations this formula is rarely used due to the need to evaluate derivatives 
of arbitrary high orders. For our purpose, however, it has the advantage that it shows that a 
real Laplace transform u)(£, t) leads to a real smearing function uj(v, t a b). Moreover, one does 
not need to know the pole structure of cl>(£,t) in the complex £-plane, which is in general 
not known. 

The result ( 1241) can be mildly generalized to the smearing of composite Hamiltonians 
vH(p,x) = vHi{p,x) + H 2 (p,x) as long as [i^i,i7 2 ] = 0. At the same time it should be 
stressed that the entire approach fails for time-dependent Hamiltonians. 

To end this section we note that if the smearing distribution uj(v,t) would have included 
negative v- values, the integral in Eq. (fl8"|) would have been replaced with f Q \— > f^ 00 - Then 
a two-sided Laplace transformation would have brought us to the same equations (j25]) - (JM|) 
as before, and the associated general u(v,t) would been recovered via the inverse of the 
two-sided Laplace transform. By restricting ourselves to v > 0, all calculations become 
simpler. 
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IV. EXPLICIT REPRESENTATION OF H(p, x) 



We now determine the Hamiltonian H(p, x) explicitly in the general case. For this we 
must first make sure that the path integral for the initial distribution in 

P v (x b ,t b \x a ,t a ) = f X(tb) Xb VxVpe^ d ^ p± - vH ^\ (30) 

Jx(t a )=x a 

is properly defined. The classical Hamiltonian H(p, x) must be set up in such a way that 
the Hamiltonian operator driving the time evolution 



d t P v (x,t\x a ,ta) = -vH{p,x)P v (x,t\x a ,t a ) 



(31) 



has all momentum operators p = —id x to the left of the x-variables. Only then can one 
guarantee the probability conservation law for P(x b , h\x a , t a ): 

d.4P.(*,H*.,tJ = - J dzvH(p, X )P„( X ,t M . (32) 

The right-hand side vanishes after a partial integration. The relation between H(p, x) and 
H(p,x) is explained in Ref. lj. 

With the help of Post's inversion formula ( |29i) we now rewrite the smeared conditional 
probability (j4]) in the form 

P(x b7 t b ;x a ,t a ) = lim^—L- d£ e^dl^Uba) 
fc-»oo [k - 1)! J Q ? 

rx(t b )=x b 

x / VxVp e L a dT(ipx-kH/o^ (33) 

J x(ta)=Xa 

Inserting here Eq. (fT9l . the integration over £ turns into 

dvv k u(v,t ba ) d< < ; '' ( r - \ VxVp e^ dT [ip± - kH/ t ] . (34) 

JO Jx{t a )=x a 

The path integral can be written in Dirac operator form (fTOj) as 



{x b \e 



-ktba H/S 



(35) 



so that fl34l) becomes 



dv v k u(v,t ba ) (x b 



1' a 



(36) 
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The ^-integral yields a Bessel function of the second type 

2 (kt ba H/y) ^ K k (2 \Jkyt ba H^j . (37) 
Applying the limiting form of the Bessel function for large index 

K k (Vkx) ~ l - k - k / 2 Y(k)(^Ye- x2 l\ (38) 

we can cast (|33|) to 



P(x b ,t b \x a ,t a ) = (x b 



00 



Xa). (39) 



dvu(v,t ba )e~ vtb « H 
We now use Eqs. ( IT9|) and (124)) to rewrite this as 

(x fe p(t feaJ ff,t ba )|x a ) = (x&l^if, l) w |x a ). (40) 
With Eq. (J2ZD, this becomes 

P(x 6 ,t 6 |x a ,t a ) = (xble-^^lxa), (41) 
which can be expressed as a path integral 

P(x b ,t b \x a ,t a ) = I ^ VxVp e^ dt[ip± - F ^ H)] . (42) 

J x(ta)=X a 

Here F C \(H) denotes the classical function of the energy which makes the path integral ( )42l) 
equal to the operator expression F(H) in fj4Tl) . The construction of fj42l) is highly nontrivial 
since one must ensure that the path integral leads to the correct operator order in u(H, l) tba . 
The task is simple only if H depends only on p. Then the ordering problem disappears and 
F cl (H(p)) = F{H(p)). 

There are other instances in which the operator ordering can be uniquely assigned. In 
Appendix B we discuss one such example. 

A few observations are useful concerning he nature of F(H). First, the normalization 
condition 

POO 

1=1 dv u(v,t) = u(0,t) , (43) 
Jo 

implies that F(0) = 0. Second, the necessary positivity of w(v, t) implies, via Post's formula, 
that all derivatives {—l) k d k uj{^ ) t)/d^ k are positive for large £ and k. Taking into account 



that u(x, t) > we conclude that u)(£, t) must be decreasing and convex for large £ and any t. 
Asymptotic decrease and convexity of e~ F ^ x ^^ ensure that F(H) must be a monotonically 
increasing function of H for large spectral values of H. Third, the real nature of u(x,t) 
makes F(H) a real function of H. This is in contrast to quantum-mechanical path integrals 
where the smearing function u(v,t) is not necessarily real. 



V. KRAMERS-MOYAL EXPANSIONS 



Let us study the implications of the above smearing procedure upon the time evolution 



lie 



equations for the conditional probabilities P v (x b , t b \x a , t a ) and P(x b , h\x a , t a ) which have t 
general form (l3Tj) . In statistical physics these are called Kramers- Moyal (K-M) equations 1 
6_|, 23]. The negative time evolution operator —vH(p,x) is called Kramers-Moyal operators 
h v (—d x ,x). Thus Eq. (IHTj) for P v (x b , t b \x a , t a ) and an analogous equation for P(x b , t b \x a , t a ) 
are written as 

d tb P v (x b ,t b \x a ,t a ) = h v P v (x b ,t b \x a ,t a ) , (44) 
d tb P(x b ,t b \x a ,t a ) = L P(x b ,t b \x a ,t a ) . (45) 

The Kramers-Moyal operator L„ has the expansion 

oo 

h v (-d Xb ,x b ) = i-dxj Di n \x b ,t b ) , (46) 

n=l 

whose coefficients Di n \x,t) are equal to the moments of the short-time transition probabil- 
ities: 

1 1 f°° 

D<?\x,t) = -Km-\ dy(y-x) n P v (y,t + r\x,t). (47) 

Inserting Eq. (13T|) . these can also be calculated from the formula 

1 1 /*°° 

Di n \x,t) = -\im- dy(y-x) n (y\e- vHT \x). (48) 

Til r^O T 

The same equation holds for P(x b , t b \x a , t a ) with the replacement vH — > H, h v — > L, 
and D^ n \x,t) -> SW(i,t). 

Our smearing procedure can be recast in the language of K-M equations. We show that 
the two equations (jlij) and (j4"5l) can be replaced by an equivalent pair of K-M equations. 
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First we observe that Eq. (IT7|) can be rewritten as 

POO 

u{z,t)= dz 1 P ul (z,t\z',t l )uj(z l 1 t l ), (49) 
Jo 



with the conditional probability 



P u (z, t\z f , t>) = -±^6{tz-t'z') u> i^f, t-?j, (50) 
which satisfies the initial condition 

limP UJ (z,t + T\z',t) = 5 + ( y z-z'). (51) 



Moreover, expressing u(z,t) in Eq. ( 1491 in terms of Pu(x b , t b \x a , t a ) via ( 1501) . we find that 
-Fw^b) ^fol^a, t a ) satisfies a C-K equation 

Puj(x b ,t b \x a ,t a ) = J dx P U) (x b ,t b \x,t)P u (x,t\x a ,t a ). (52) 

This implies that Pu(x b , t b \x a , t a ) obeys a K-M time evolution equation 

d tba Pu(v b ,t b \v a ,t a ) = LuPu(v b ,t b \v a ,t a ) , (53) 

with the K-M operator 

oo 

L^^KfxW(,,g. (54) 

n=l 

The expansion coefficients are obtained by analogy with f|4Tj) from the short-time limits 

i r°° 

K^ n) (v,t) = lim — / di/(v / --u) n P (i ,(« / ,t + rkt) 
r^o n!r J_ oc 



' — oo 

n poo 



= i ir[1 ^_ ( — I — ) / dv'(v'-v) n uj(v',T). (55) 
r-^onlr \t ba + rj 7 

A crucial observation for the further development is that due to the equality of the 
Kernel in the time evolution equation ( 1531) and ( 1491) . the same K-M operator governs the 
time evolution of u>(v, t ba ): 

d tba u(v, tba) = 1^ t^) , (56) 

In Appendix A we show that this equation can also be derived directly from the K-M 
equations (jSj) and f|4"5j) . 
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The smearing procedure of the path integral is equivalent to replacing the pair of K-M 
equations (pH|) and fj45|) by the pair (jHJ) and fl56|) . This separates the dynamics of the 
smearing distribution from the dynamics of the transition amplitude. This may serve as a 



24j . in superstatistics [if]], or in 



convenient starting point, for instance in quantum optics 
numerous option pricing models (see, e.g., Ref. jl] and citations therein). 

For applications it is useful to remember Pawula's theorem 25j, according to which the 
coefficients of the expansions of K-M operators are either all nonzero or, if there is a finite 
number of them, they can be nonzero only up to n = 2. This follows from the necessary 
positivity of the probabilities. If one artificially truncates Di n \x,t) in ( H6|) or K^ n \x,t) at 
some n ^ 3, then the ensuing transition probabilities always develop negative values, at 
least for sufficiently short times. This is the basic reason why phenomenological models for 
K-M operators go usually only up to n = 2. 

Consider such a truncated model. Then the K-M equations (jHJ) and (156]) reduce to the 
Fokker-Planck equations 

dtu(v,t) = h F fu;(v,t), (57) 
d tb P v (x b ,t b \x a ,t a ) = h FP P v (x b ,t b \x a ,t a ) , (58) 

with 

I/J = -d v K^(v,t) + d 2 v K^(v,t), (59) 
Lf = ~d Xh ( Xb7 t b ) + dl b (x b , t b ). (60) 
Furthermore, Eqs. (157|) and fl58|) allow us to find two coupled stochastic processes described 



by the two coupled Ito stochastic differential equations 



dx b = Di 1 \x b ,t b )dt b + \/2Di J 2 \x b ,t b )dW l , (61) 



d^ = K^\v, t ba ) dt ba + ^2K&(v, t ba ) dW 2 . (62) 
Here Wi(t b ) and W2(t ba ) are Wiener processes, i.e., Gaussian random walks. 



VI. SIMPLE EXAMPLES 



To demonstrate usefulness of the superposition procedure we now discuss some important 
classes of G(x)-functions in Eq. (1241) . 
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i.) We start with the trivial choice 



G{x) = e~ ax+b , beR; ae 



^0 ' 



which gives 



and consequently 



u((,t) = e 



-a(+bt 



uj(v,t) = e bt S(v-a) 



(63) 



(64) 



(65) 



By requiring that uj is normalized to 1 for any t we have 6 = 0, i.e., no smearing distribution. 
In this case the Hamiltonian H(p, x) = aH(p, x). Note that a is basically the averaged value 
of v over the 5-function distribution. 



ii.) A less trivial choice of G(x) is 



G(x) 



a G R + ; b, c G 



which gives 



x + b 



w(C,*) 



^0 ) 



at 



C + bt 



ct 



leading thus to 



uj(v, t) 



{at) ct e- btv v ct - 1 



T(ct) 



(66) 



(67) 



(68) 



Further restriction on the coefficients is obtained by requiring normalizability of u(v, t). The 
normalization condition F(0) = can be fulfilled in two ways. Either we set c = 0, in which 
case u(v, t) = 5(v) and H(p, x) = 0, or we assume that a = b. In the latter case 



u>(v, t) 



(bt) ct v 



ct„,ct—l 



-btv 



tance in financial data analysis [1| 



(69) 



net) ' 

i.e., it corresponds to the Gamma distribution [4j, [lOj fu,ct{ v ) which is of a particular impor- 



26] and superstatistics [27J. In the special case when 



c = db = da and b — > oo we obtain that u(v,t) = 5(v — d). The Hamiltonian H associated 
with the distribution (1591 reads 



H(p, x) = vb 



logl'^+l 



(70) 
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where v = c/b is the mean of u(v,t). In particular, for H = p 2 /2 we obtain the smearing 
relation between path integrals 

VxVp e / to b dr^-^log(pV26 + l)] 

% (j^a )— &a 

poo rx(t b )=x b t 

= dvu(v,t ba ) VxVpe^ dT{ip± - vp/2 \ (71) 

JO Jx(t a )=x a 

Inserting for the path integral on the right-hand side the result ([151) . the integral over v 
reads explicitly 

"°° (bt ba ) ctba v ctba - 1 I 1 



and yields 



dv K ba \ J^—e- bt »° v e-</ 2vt *«, (72) 

T{ct ba ) \l 2nt ba v K ' 



K ^^vk^(^f) ' (73) 



ia, ii6i, 12 



In Fourier space, the superposition (ITT]) leads to a Tsallis distribution 

Another interesting consequence arises when we consider the ensuing Ito stochastic 
equations. To this end we use the Hamiltonian from Refs. [l[ [lo| which has the form 
p 2 /2 + ip(r/v — 1/2), with r being a constant. The corresponding drift and diffusion coeffi- 
cients Dy and are easily calculable giving 

D$-\x,t b ) = (r- V -), D?(x,t b ) = \. (74) 

The coefficients are 

K W (v,t ba ) = ^(- h ~v) = ^-(v-v), 
t ba \o / t ba 

K^(v,t ba ) = —4-, n>2 - (75) 
v ; t% a nb n K J 

Hence we find the two coupled Ito equations (IBTj) and (|62|) in the form 

dx b = {r-D^b +VvdW 1 , 

dv = ^-(v-v) dt ba + -L/f dW 2 . (76) 

ha Ha V 

One may now view x b as a logarithm of a stock price S, and t> and r as the corresponding 
variance and drift. If, in addition, we replace for large t ab the quantity y/ti with y/v, the 
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systems ( 1751) reduces to 



dS = rSdtb +y/vSdWi 



dv 



v) dt ba + e^/v dW 2 



(77) 



where 7 = l/t& a and £ = l/(6t& a ). The system of equations (177)) constitute Heston's stochas- 
tic volatility model [261]. The parameters u, 7 and e can be then interpreted as the long-time 
average, the drift of the variance and the volatility of the variance, respectively. Heston's 
model may be used in quantitative finance to evaluate, for instance, the price of options (for 
further reference on this model see, e.g., Ref. [l] and citations therein). 

iii.) As a third example we consider 



G(x) 



This implies 



and 



u>((,t) 



u(v, t) 



a G 



ae *> 



(78) 



(79) 



(80) 



2V^V f 

Image function uj((,t) already fulfils the normalized condition F(Q) = 0. In the literature, 
( l80l) is known as the Weibull distribution of order 1 (e.g., Refs.j^, 32]) u(v, t) = w(v; 1, a 2 t/2). 
The smearing distribution (1801) has also important applications in the so called inverse 



X superstatistics : 



Because integral over v a u(v,t) does not exist for a > 1/2, the 



distribution (IHUl) does not have moments. The Hamiltonian H(p, x) = a[\J H(p, x)] c i. In 



particular, this implies an identity 



£\ rx(t b )=x 



dvwlv; 1, - 



x(t 6 )=x 6 



(81) 

'x(t a )=X a 

The right-hand side represents the path integral for free relativistic particle in the Newton- 
Wigner representation [30|, l31|. Previously, this has been evaluated by group path integra- 
tion [33| . With our smearing method, we can obtain the same result much faster by a direct 
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calculation of the left-hand-side of (IHTj) . Due the the quadratic nature of the Hamiltonian 
we obtain immediately in D space dimensions: 

fx(t () )=x i , 



r X 'l)xDp e ^ d K ip ^ c v / P i +™ i ^) 

Jx.(t„,)=Xa 



'x(i a )=x a 

— „ „ x 2 / 1 X-D/2 
e *■» -vm 2 cH ba S 2 — ' I 

dv =e 4vc2t ba 



2 H£ 



4vc 2 irt ba 



D+l 

with 7 = (1 + xj^/c 2 ^) -1 ' 2 , Xfe a = xt — x a . The result in Eq. ( 1821) agrees, of course, with 
those of earlier authors [33[ ] . 

If Eq. ( 182]) is taken to the limit m = 0, it reads 

M>p e/- 6 IpD = at r (2±1) . (83) 

This is of a particular importance in econophysics as it can be directly related to a Levy 
noise distribution of order 1, see, e.g. Ref. [1]. 



VII. CONCLUSIONS AND OUTLOOKS 



In this paper we have identified the most general class of continuous smearing distributions 
lu(v, t) defined on M + x M + for which superpositions of Marcovian processes remain Marko- 
vian. This insight was used to rephrase the original problem of computing a complicated 
conditional probability P(xb,tb\x a ,t a ) as a problem solving two coupled Kramer s-Moyal 
equations, one for a simple basic distribution, and the other for the smearing distribution 
u(v,t). If the associated Kramers-Moyal expansions are truncate after n = 2, they become 
Fokker-Planck equations whose respective sample paths follow two coupled stochastic equa- 
tions a la Ito. As an application of this decomposition procedure we have demonstrated 
how the log-normal fluctuations of a stock price with a Gamma smearing distribution of the 
variances lead to the celebrated Heston model of stochastic volatility. 

As a second application we have shown that the ensuing relation between smeared path 
integrals and non-smeared ones permits one to solve exactly (and often quite fast) relatively 
large classes of path integrals. In this connection we have discussed in some detail two 
important situations: path integrals arising in the framework of generalized statistics of 
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Tsallis 



28] and a world-line representation of the euclidean Newton- Wigner propagator. 



The latter situation can easily accommodate relativistic systems with gauge potentials via 
minimal coupling. Useful applications are expected in the field of world-line representations 
of effective actions Q, Q . 

For simplicity we have considered smearing distributions for path integrals with only 
bosonic degrees of freedom. There is no problem in extending our procedure to path integrals 
with Grassmann (i.e., fermionic) variables, and to more general initial functional integrals. 
Such extensions may be useful in polymer physics, in particular in theory of self-avoiding 
chains, and in their field-theoretic treatments [ll. 136]. 

Our procedure can also be applied to quantum mechanics with only small modifications. 
There the C-K relation is the composition law for time translation amplitudes reflecting the 
semigroup property of the evolution operator e - l f dtH . The smearing functions acts then 
upon probability amplitudes rather than conditional probabilities, in which case they may 
be complex rather than positive. 

When the manuscript was finished we have became aware of Refs. {itI . I^ . devoted to 
superstatistics. There exists a conceptual overlap between the present paper and the above 
mentioned ones. In fact, practitioners in superstatistics can quickly benefit from our paper 
by substituting for our v their intensive parameter (3 and for our smearing distribution 
u(v,t) their probability density /(/?). In superstatistical setting one could also call our 
smearing procedure as "the superstatistical average" (though our distributions are generally 
time dependent). In any case, it seems to us that important things are worth repeating 
several times using different framework and different words. For a more detailed discussion 
of the superstatistics paradigm, the reader is referred to Refs. 

HH, El- 
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Appendix A 



It is instructive to check that Eq. (I58|) can be derived directly from the original K-M 
equations (jSJ) and (HS]) for P v (x b ,t b \x a ,t a ) and P(x b ,t b \x a ,t a ). This provides an important 
cross-check for the complete equivalence between both systems of K-M equations. To this 
end, we multiply (jHJ) by uj{y,t ba ) and integrate over v. This leads to 

POO 

d tb P(x b ,t b \x a ,t a ) - / dvP v (x b} t b \x a} t a )d tba uj(v,t ba ) 







poo 

/ dvu(v,t ba )'L v P v (x b ,t b \x a ,t a ). (84) 
Jo 



The first line has been rewritten using the chain rule for derivatives, while in the second line 
we have used the fact that P v (x b ,t b \x a ,t a ) fulfills the K-M equation fjl4|) . We now insert in 
the second line a trivial unit integral J °° dv' u(v',r) = 1, where r is a positive infinitesimal 
time. Then we subtract and add a term L P(x b , t b \x a , t a ), i.e., the right-hand side of the 
K-M equation fT4"5l) . Thus the second line in (1511) becomes 



OO POO 



dv / dv'u(v,tba)u}(v' ,t)[L v -L, v >}P v (x b ,t b \x a ,t a ) 
o Jo 

+ LP(x b ,t b \x a ,t a ). (85) 

Eq. (|85p is a convenient short-hand version of more involved form where the r — > limit 
acts both on uj(v, t) and L„. For instance, 

POO 

. . . dv'u(v', t)Iv . . . , (86) 
Jo 

actually means 

°° 1 POO POO 

...V—/ dv'u{v',r){-d Xb ) k dy{y-x b ) k P v ,{y,t b + T,x b ,t b ) ... , (87) 

k=l Tl£ -Jo J-oo 

and similarly for the L„-term. With this in mind we now Taylor-expand L„/ around v and 
obtain 

-V—/ dw / dv'u(v,t ba )uj(v',T)(v' - v) n [d™L v ) P v (x b ,t b \x a ,t a ) 

n=l n! 

+ LP(x b ,t fe |x a ,t a ). (88) 

To compute the derivatives d™h v we use the representation for a more general unsmeared 
Hamiltonian operator H(p, x): 

P v (x b ,t b \x a ,t a ) = e- vt »« A ^5(x b -x a ) = e- vtb ^^ x ^5(x b -x a ), (89) 
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from which directly follows that 



d2P v (y,t b + T\x b ,t b ) = (-rrH n (p,y)e- VT ^5(y-x b ) 

) n [(#) n (p b ,x b )]*e- VT &fo*W6{y-x b ). (90) 



-T 



Appearance of the term [W(p,x)]* in the second lines is a consequence of the identity 
(x b \A(p,x)\x a ) = (x a \A*(p, x)\x b )* which is valid for any operator A. We can now simplify 
(188)) with some algebra involving product rules of derivatives and 5-functions. Abbreviating 
H(p,x) by H x , this gives 

dv dv'u(v,t ba )tu(v',T)(v' -v) n d:P v (x b ,t b + T\x a ,t a ) 

nh \tba + rj J J 

+ LP(x b ,t b \x a ,t a ). (91) 
In deriving the latter we have utilized the identity 

Jdy(y-x b ) k [(Hl) n Ye~ VT 5{y - x b ) = J ' dy5(y - x b )H^-^(y - x b ) k , (92) 



and the fact that 



V 



Hu e -rvH y = ( __L ) e~ TvH y . (93) 



Extending for a moment the sum over n in (leg)) to start from n — 0, the relevant terms 
entering the calculation of [<9™L tI ] P v (x b , t b \x a , t a ) in Eq. (|9T)) have the form 

E k\ 5{V ~ Xb >~ TVHv ( y ~ ^) k e- tb ^5(x b - x a ) 

k=0 

= E lz # L 6 (y - ^)h/(- <T ) - y\ k e- tbaVHy Kv - x a ) 

k=0 

= e [v - v{ - iT ^5{y-x b )e- tbaVH »5{y-x a ). (94) 

We now apply the shift operator e^ y ~ y (~ lT ^ dx b to the ^-variable in the 5-function and cast 
the resulting S(x b — y(—ir)) into Dirac's bra-ket form 

6(x b - y{-ir)) = (x b \y(-ir)) = (x b \e~™ 6 \y) . (95) 
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If we now perform the integral over y in (|92|) and use the fact that 

C^-J (x b \e-^ +t ^\x a ) = (r^hf) (* b \e- ViT+tba) %a) , 

we obtain the announced result (I9~T]) . We now use the definition of the moments in the 
second line of Eq. (1551) . and the fact that the surface terms in the partial integrations over 
v' vanish due to the positivity of the real part of the Hamiltonian spectrum, to rewrite the 
first term in ( 19T1) as 

OO /»00 

-V / dv(-d v ) n [K (n \v } t ba )uj(v,t ba )]P v (x b ,t b \x a} t a ). 

n=l J( > 

Remembering Eq. (1541) . this can be written as 

-J dvP v (x b ,t i \xa,t a )K>u(v,t ba ). (96) 
If we now take into account the evolution equation (PHI) , we obtain from Eq. (1M1) the equation 

/oo /* 
dv P v (x b ,t b \x a ,t a ) d tba u{v,t ba ) = J dv Pvix^tblx^Qh^uiv^ha). (97) 

This result must hold for any distribution P v (x b ,t b \x a ,t a ), thus proving the K-M equation 
(PI) for u(v,t). 



Appendix B 

Here we show how to define uniquely the classical Hamiltonian F c \(H(p, x)) if the initial 
Hamiltonian H(p, x) depends not only on p but also on x, This will be done with the help 
of the K-M equations. 

We start with an observation that when the unsmeared K-M operator has no explicit 
time dependence, i.e., when 

oo 

K(-d x ,x) = ^(-d x ) n Di n \x), (98) 

n=l 

then the unsmeared classical Hamiltonian H c \(p,x) has the form 

oo 

H cl (p,x) = -J2(-ip) n Di n \x), (99) 

n=l 
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provided we define the path integral by time slicing in the post-point form. This 
is done by rewriting the short-time matrix elements (x n \ exp(— Hr)\x n -i) as integrals 
f(dp n /2ir)(x n \p n ){p n \ exp(— Hr)\x n -\) and defining the classical Hamiltonian via the iden- 
tity 

(p n |exp(- J ffr)|x n _i) = e-rH^Pn^) e -i PnXn . K 

Then to order 0(t 2 ) the Hamiltonian H c \(p,x) coincides with ( 1991) . Full discussion of this 
relation can be found in Ref. The K-M equation is basically a Schrodinger-type equation 
with the non-hermitian Hamiltonian H(p,x) = 'Yl in {—ip) n D^ l \x) = — L„(— d x , x). 

Let us now see how the K-M equation looks for the smeared distribution. This will allow 
us to identify the smeared classical Hamiltonian F C \(H). To this end we note that from the 
definition fH7|) the smeared K-M coefficients can be written as 

1 i /*oo /»oo 

D {n \x) = — lim-y dvu(v,r) J dy (y - x) n (y\e- vflr \x) 

= ^lim^ J dy{y-x) n (y\e- F(H)r \x). (101) 

Furthermore, we can insert in front of e~ F ^ T a completeness relation J dp \p)(p\ = 1, and 
make use of the fact that only terms up to order 0(t) are relevant in e~ F ^ T . This leads 
directly to 

i-lim- / dpdy (y-x) n (y\p)(p\e~ F{fl)T \x) = —lira- f y n e ipy e~ F ^ H)T . (102) 

Here F C \(H) is obtained from F(H) by using the commutator [p,x] = —i to move all p's 
to the left of all x's and dropping the hats. Eq. (11021) can be further simplified if we write 
F C \(H) = Y^nP n fn( x ) an d the ^/-integral as the n-th derivative of 5-function, i.e. 



D {n \x) = ^-Mim- / dp e - rF ^ H) d?8(p) 
i n 1 

= 9 n P [e- TF « (H) ] U = (-i) n fn(x) . (103) 
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